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An accretion flow is necessarily transonic around a black hole. However, 
around a neutron star it may or may not be transonic, depending on the inner 
disk boundary conditions influenced by the neutron star. I wiU discuss various 
transonic behavior of the disk fluid in general relativistic (or pseudo general 
relativistic) framework. I will address that there are four types of sonic/critical 
point possible to form in an accretion disk. It will be shown that how the fluid 
properties including location of sonic points vary with angular momentum of 
the compact object which controls the overall disk dynamics and outflows. 

1 Introduction 

Accretion disks serve to identify compact objects, mostly black holes and neu- 
tron stars, in the universe. The most common way to understand its formation 
is in a binary system where matter is pulled off a companion star and settles 
on to the compact object in the form of a disk. As practically black holes 
and neutron stars cannot be seen, by detecting and analyzing light rays off 
an accretion disk one can understand the properties of the central compact 
object. Other examples of the accretion disk are the protoplanetary disks, 
disks around active galactic nuclei, in star-forming systems, and in quiescent 
cataclysmic variables etc. 

The molecular viscosity in an accretion disk is rather low. Shakura & Sun- 
yaev [T] proposed that there is a significant turbulent viscosity explaining 
transport of matter inwards and angular momentum outwards with the Ke- 
plerian angular momentum profile 2 . However, until the work by Balbus & 
Hawley [3^ origin of the instability and plausible turbulence was not under- 
stood. Applying the idea of Magneto-Rotational-Instability (MRI) described 
by Velikhov ^ and Chandrasekhar 'W , they showed that the Keplerian disk 
flow may exhibit unstable modes under perturbation in the presence of a weak 
magnetic field. This plausibly generates turbulence which results transport. 
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However, there are several accretion disk systems where gas is electricahy neu- 
tral in charge and thus the magnetic field does not couple with the system. 
Examples of such system are protoplanetary and star-forming disks, disks 
around active galactic nuclei and in quiescent cataclysmic variables. Hence 
in these systems MRI is expected not to work to generate turbulence. While 
the origin of the transport in such systems is still an ill understood problem, 
some mechanisms have been proposed by various groups including the present 
author [6l [H [HI [9l [IHl El [12] . 

However, the problems of the transport and then the viscosity are not se- 
vere if the disk angular momentum profile deviates from the Keplerian type to 
sub-Keplerian type. In a sub-Keplerian disk, the gravitational force dominates 
over the centrifugal force resulting in a strong advective component of the in- 
falling matter [13l [lU [15] . Due to the dominance of the gravity, the matter 
even with a constant angular momentum may transport inwards easily. An 
accretion disk with a strong advective component, namely the advective accre- 
tion disk, is hotter than the Keplerian disk which explains successfully several 
systems e.g. Sgr A*, GRS 1515-fl05 which are expected to be consisting of 
hot disks. Therefore, such an accreting system may not be of pure Keplerian 
type. At least close to the compact object it is an advective accretion disk. In 
a sub-Keplerian disk, at far away off the black hole, the speed of the infalling 
matter is close to zero when it is practically out of the black hole's influence. 
However, temperature of the system is finite resulting in finite sound speed at 
that radius. On the other hand, matter speed at the black hole horizon reaches 
the speed of light (c) , while, by causality condition, sound speed can be at the 
most c/-\/3- Therefore, there must be a location where matter speed crosses 
local sound speed making the flow transonic. However, the second condition 
namely the inner boundary condition not necessarily be satisfled in the flow 
around a neutron star. Therefore an accretion flow around a neutron star may 
or may not be transonic. 

In the present article, I mainly concentrate on the sub-Keplerian accre- 
tion disk with possible multitransonic flow. If the flow is transonic, then the 
disk dynamics and corresponding outflow are influenced by the sonic loca- 
tions/points [m [ini [m [m [I3. Depending upon the flow energy (entropy) 
and angular momentum, a disk may exhibit one to four sonic points, through 
which all however, matter may not pass. Sonic points in accretion disks play 
role similar to the critical points in simple harmonic oscillators, particularly 
to control the dynamics of the system. Indeed an accretion disk can be looked 
upon as a damped harmonic oscillator [20j . 

I organize the paper as follows. In the next section, I discuss accretion 
flow, flrst the spherical one and then the disk flow, and formation of possible 
sonic points therein in the absence of any energy dissipation. In §3, I analyze 
how the rotation of the compact object affects the sonic points. Subsequently, 
for completeness, I present the generalized set of accretion disk equations 
including possible dissipation effects in §4, without going into their solutions. 
Finally, I summarize in §5. 
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2 Accretion flow and formation of sonic points 

First I discuss simple spherical accretion namely Bondi flow. Subsequently, I 
include angular momentum into the equations set describing disk accretion. I 
show that while the Bondi flow around a nonrotating compact star exhibits 
single sonic point the disk flow may exhibit multiple sonic points. 

As gravity is expected to be very strong close to the compact object, in 
principle I should describe the flow system by the set of general relativistic 
equations. For a nonrotating compact object equations should be written in 
the Schwarzschild geometry and for a rotating compact object they should be 
in the Kerr geometry. However, that might hide transparency of the descrip- 
tion, as the full general relativistic set of equations is so cumbersome that it 
is difhcult to relate the terms with the physics they carry. Indeed under some 
occasions a full general relativistic description is not required. Therefore, I de- 
scribe the system by pseudo-Newtonian approach. Here one uses Newtonian 
set of equations only along with a modified gravitational force/potential such 
that it mimics the general relativistic features approximately to describe the 
system. 



2.1 Bondi flow 



This happens for an isolated star when matter falls onto it from all directions, 
resulting in a spherical accretion. Therefore, to describe the system, we con- 
sider the spherical polar coordinate system where only nonzero component of 
velocity is Vr = v. Hence the equations describing the steady-state flow with 
negligible viscosity are given as: 

dv IdP 



~^"7" (pvx^) = 0, hence pvx^ — constant = Mac, (2) 
x-^ dx 

where all the variables are expressed in dimcnsionless units, v is velocity in 
unit of light speed, x = r/rg where r is radial coordinate and Vg = GM^jc? ^ 
Mg,c,G are respectively mass of the central star, speed of light, Newton's 
gravitation constant, P, p are respectively corresponding dimcnsionless pres- 
sure, density of the flow, Mac is accretion rate and Fg is gravitational force 
given by 

Fg — — : Newtonian, (3) 
x^ 

pseudo — Newtonian, Schwarzschild geometry pT| , (4) 



(a; -2)2 

{x^ - 2j^ + i'^f 
a;3(Vi(x-2)+j)2 



pseudo — Newtonian, Kerr geometry [22] , (5) 
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where j is specific angular momentum of tlie compact object. Integrating 
equation ([T]) for an adiabatic flow to a nonrotating compact object I obtain 

E=^v'+na'-^, (6) 

wliere polytropic index n = 1/(7 — 1) and 7 is defined as 7 = a^p/P with 
P = Kp"', K is gas constant and a is sound speed. If the Mach number of 
the flow is defined as M = v/a, then at a constant energy E, it is clearly 
understood from ^ that M — x trajectory is always hyperbolic-type [53] 
with a sonic point at M = 1 shown in Fig. [1] ^20]. The solution marked 'A' 
by an arrow indicates accretion and marked 'C" as wind. None of the other 
solutions are physical to describe accretion and wind. However, for the flow 
around a rotating compact object, there might be more than one sonic point 
but I will not discuss them here, this can be checked with a pseudo-Newtonian 
potential proposed for the Kerr geometry [52]. 



2.2 Disk accretion flow 



Now I consider a binary system when the compact object is closely associated 
with its binary companion star and pulling matter off it. The gravitational 
force acts on the donor star radially. On the other hand, the star is rotating 
and hence the matter detached off it due to gravitational pull has angular 
momentum. As a result the matter infalls towards the compact object in a 
spiral path forming a disk around it called the accretion disk. Therefore, the 
radial momentum balance equation describing disk dynamics in steady-state 
in the absence of any energy dissipation is [I8j 

dv IdP X"^ 

"^Tx^-pl^-V^^"^^-'^ 

where A is specific angular momentum of the infall which is a conserved quan- 
tity in a nondissipative fiow. The corresponding vertically averaged mass con- 
servation reads as 

— ^ (xSv) — 0, hence — 4,'n:xEv = constant = Mac, (8) 
X dx 

where S ~ h{x)p is the (vertically integrated) column density [24j with p 
is density at the equatorial plane, h{x) = a^Jx/Fg is the half-thickness of 
the disk computed from the assumption of vertical equilibrium. Now integrat- 
ing ((Tj) for an adiabatic fiow towards a nonrotating compact object I obtain 
specific energy of the flow 
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Fig. 1. Upper panel: Bondi flow: Mach number as a function of logarithmic radial 
coordinate for various energy of the flow. The point 'B' is sonic point and the solution 
marked by the arrow 'A' indicates accretion and 'C'' indicates wind. Lower panel: 
Disk flow: Same as Bondi flow; [20] . 
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which is a conserved quantity for the nondissipative system. It is seen from 
^ that at far away from the black hole and close to its event horizon lo- 
cated at a; = 2, last term (gravitational potential energy) on the right hand 
side dominates over the last but one term (centrifugal energy) resulting in a 
hyperbolic-type M — x trajectory as shown in Fig. [1] [16l [20]. On the other 
hand, at an intermediate location, with an appropriate value of A, centrifu- 
gal energy may dominate over gravitational potential energy resulting in an 
ellipse-type M — x trajectory as shown in Fig. [1] 



2.3 Formation of sonic points in disk accretion 



Combining ([T]) and ([8]) I obtain 



dv 
dx 



F, 



^ 7+1 \ x Fg dx ) 



(7 + 1)1' 



N 
D' 



(10) 



At the sonic point x ~ Xc, D = and thus to have continuous dv/dx at the 
same location N must vanish. Hence, applying N = D = a.t x = Xc to PH)) . 
I obtain velocity, sound speed, specific energy at Xc 



7 + 1' 



hence Mach number Mr 



7 + 1' 



(7 + l)(-J-^^.c 



27 



7-1 



A2 



1 fdF„ 



_3_ 

Xr 



J_ _ J. 

I also compute a quantity that carries information of entropy 



(11) 



(12) 



which is a useful quantity to study the sonic point behavior. Now to obtain 
solutions for the accretion disk one has to integrate PH]) with a proper bound- 
ary condition. As our interest is in transonic solutions, the matter necessarily 
passes through the sonic location and this I consider as our staring point of the 
integration. As the system under consideration is in the steady-state, integrat- 
ing from the sonic point to outwards (inwards) and vice versa are equivalent. 
However, dv/dx given in PT!]) is of 0/0 form. Therefore, applying I'Hospital's 
rule, I obtain dv/dx at x = Xc 



dx)^ 2 A 



(13) 



where. 
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A 
B 

C 



= 1 



(7 + IK' vlin + W 



Veil + If 
^2 



3 



dx 



7 + 
3A2 



1 



dx J 



dFg 

dx 



1 fd^F„ 



F, 



«g(7-l) 
(7+1)' 



dx'^ 



Fr„ 



dx 



(14) 



Now the value of discriminant T) = — AAC determines the type of sonic 
point. Clearly T) depends on the sonic location which can be computed for a 
given Ec from (fTTj) . In fact, the expression for Ec can be rewritten as a fourth 
order polynomial of Xc. Hence, Ec and A are the input parameters determining 
flow behavior for a given system. 

There are four different types of sonic point shown in Fig. [2] pS] l20j classi- 
fied according to the trajectory around it. (1) When I? < with S = 0, sonic 
point is center-type (elliptical trajectories in Fig. [T]). (2) I? < with B ^ 
gives spiral-type sonic point. For T) > with (3) AC > it is nodal-type, and 
(4) AC < 0, saddle-type (hyperbolic trajectories in Figs.[T]). When AC — for 
T) > 0, sonic points are called straight line. 




C«ntnr SvfnsT Sod^te No^pl 



Fig. 2. Trajectories around various critical points where abscissa and ordinate are 
radial coordinate and Mach number respectively; [201 125) . 



3 Sonic point analysis in disks around rotating compact 
objects 

The energy and information of entropy at the sonic point Ec and A^c as 
functions of sonic location Xc given by (fTTj) and p2|) respectively describe the 
loci of sonic point energy and entropy of the flow. Figure [3] T8| shows their 
behavior for various choice of speciflc angular momentum of the compact 
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object. As I consider a nondissipative disk, the energy and entropy remain 
conserved throughout. The intersections of the curves by the constant energy 
and entropy hne (which are the horizontal Unes in the figure) indicate the sonic 
points of the accretion disk for that particular energy and entropy and specific 
angular momentum of the black hole (Kerr parameter). It is clearly seen that, 
at a particular energy and entropy, if the Kerr parameter increases, the region 
where sonic points form (as well as radii of marginally bound (xb) and stable 
(xs) orbits) shift (s) to a more inner region of the disk and possibility to have 
all sonic points in the disk outside event horizon increases. As an example, 
for a — 0.998, the inner edge of the accretion disk enlarges to such an extent 
that the fourth sonic point in the disk appears outside the event horizon. On 
the other hand, for retrograde orbits (counter rotating cases), all the sonic 
points come close to each other tending to overlap for a particular value of 
the Kerr parameter (when Xb and Xg move to greater radii). Therefore, as the 
value of the Kerr parameter decreases, the possibility of forming individual 
sonic points decreases and thus shrinking the region containing sonic points. 
This is very well-understood physically; as the Kerr parameter decreases, total 
angular momentum of the system decreases resulting in the disk tends to a 
Bondi-like flow, that has single saddle-type sonic point. 

On the other hand, for a particular a, if Ec and Mc decrease, then the 
possibility to have all the sonic points in the disk increases. This is understood 
from Fig. [3] that as Ec and A4c decrease, the curve is more likely to intersect 
the horizontal line. If Ec — 0.05, then number of intersection of all the curves, 
except the dashed one, by the Ec — 0.05 line is one giving only sonic point. 
The dashed curve with a — 0.998 intersects twice. However, for Ec — 0.0065 
(as considered in Fig. [3]), the dashed curve intersects four times exhibiting 
four possible sonic points in the flow and the other curves intersect thrice. For 
a detailed description, see earlier works [HI [19]. 

Note that the sonic points occurring with a negative slope of the curve 
indicate the locations of 'saddle-type' sonic point and those with a positive 
slope indicate the 'center-type' sonic point. Thus the rotation of a black hole 
plays an important role in the formation of the sonic points which are related 
to the structure of accretion disks and presumably formation of outflows and 
jets. 

4 Generalized set of equations describing a 
sub-Keplerian accretion disk 

So far, for simplicity, I have described disks without any dissipative energy. In 
principle a disk must exhibit viscous dissipation making angular momentum 
varying with disk radius. The hot disk flow with ion temperature T > lO^K 
is also expected to generate significant nuclear energy mostly via proton- 
capture reactions [26l[27j. On the other hand, significant energy is radiated 
out through the inverse-Compton, bremsstrahlung and synchrotron effects, 
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Fig. 3. Variation of energy (upper panel) and information of entropy (lower panel) 
as functions of logarithmic sonic radius when j is a parameter. Central solid curve 
is for nonrotating compact object {j — 0), while curves in regions of either side of 
it indicated by '+' and ' — ' are for prograde and and retrograde orbits respectively. 
Different curves from j = line to downwards are for j — 0.1, 0.5, 0.998 and to up- 
wards for j = —0.1, —0.5, —0.998. The horizontal line indicates the curve of constant 
energy of 0.0065 (upper panel) and of constant entropy of 5 x 10~^ (lower panel). 
Other parameters are A = 3.3, 7 — 4/3; [18] . 
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cooling the disk. In addition, significant energy is expected to be absorbed 
through endothermic nuclear reactions, mostly via dissociation of elements pSl 
\T7[ [2S] , affecting disk dynamics significantly. Although the main purpose of 
the paper is to describe the basic mechanism to form sonic points in accretion 
disks and for that the inviscid set of equations without dissipation suffices, for 
completeness here I describe the generalized set of equations including effects 
of energy dissipation. However, 1 will not go for the solutions of such equations 
set which is beyond the scope of present paper. 

Therefore, in general one should incorporate two more equations namely 
angular momentum balance and energy equation apart from ([7]) and ([8]) as 
given below. I express viscous dissipation (7+ in terms of shear stress Wx^ as 
q+ = W^^/r], where 77 is the coefficient of viscosity. Shakura & Sunyaev [1] 
parametrized shear stress in a Keplerian disk by gas pressure with a constant 
a, called Shakura-Sunyaev viscosity parameter, such that Wx^, = —aP. How- 
ever, in an advective disk, there must be a significant contribution due to ram 
pressure and thus shear stress may be read as W^^, = —a{In+iP + Inv'^ p)h{x) , 
where I consider vertically averaged disk and /„ appear due to integra- 
tion in the vertical direction |24] when P and p are pressure and density 
respectively at the equatorial plane. However, if I use this expression of shear 
stress in describing g+, it loses information of actual shear, effect due to vari- 
ation of angular velocity in the disk. The proper way of writing shear stress 
should be Wx^, — V^^i where fl is angular frequency of the flow. If I sub- 
stitute this into , then final equation to be solved will contain a nonlinear 
derivative term making it tedious to solve. Therefore, following Chakrabarti 
I express (7+ = —a{In+iP + InV^ p)h{x)x^ , where one Wx4, in the square 
is expressed by pressure and other by actual shear. Thus the angular momen- 
tum balance equation is 

dX _ I d 
dx Sx dx 
Then the energy equation may read as 

ds vh{x) (dP P dp\ + , . 

^vT— = — - Pi-— ]=Q^ -Q , (16) 

dx 1 3 — 1 \ dx p dx J 

where s is entropy density, = 9^ + <Znuc and = qic + Qhi + <Zsyn when 
Qnuc is the nuclear energy released/absorbed in the disk and Qic , 9br , (Zsyn are 
respectively energy radiated out due to inverse-Compton, bremsstrahlung, 
synchrotron effect. Following Cox & Giuli j30j, I define 

P „-, , ^ _ ^ (4-3/3)^(7-1) gas pressure 

^ 4-3/3' ' ^ (3+ 12{j~l){l- 13)''^ total pressure 

To obtain generalized solutions of an accretion disk with a significant ad- 
vective component, namely a sub-Keplerian accretion disk, one has to solve 
equations (O, ([5]), ^5)1 and simultaneously along with the set of equa- 
tions generating nuclear energy through the reactions taking place in a hot flow 
among all the isotopes therein as given in earlier works in detail [271 HH] ■ 



V— = — - [x^q^] ■ (15) 
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5 Summary 

I have described how the various types of sonic point form in an accretion 
disk. While sonic points necessarily form in disks around a black hole, it need 
not be around a neutron star. I have discussed that upto four sonic points may 
exist outside event horizon of a black hole depending on the flow parameters 
and black hole's angular momentum. However, in a zero angular momentum 
Bondi flow around a nonrotating compact object, there is only one sonic point. 
When a disk has saddle and/or nodal -type sonic points, matter passes through 
them, while it does not for center and spiral -type sonic points. Therefore, the 
last two are not physical sonic points for accretion and wind because they 
may not be the part of solutions extending from infinity to the black hole 
horizon. However, the corresponding branches may be useful to explain certain 
scenarios in disks, as argued by some authors, e.g. formation of shock helping 
to understand observed truncation of disks [3T1 [32 1 launching of outflows and 
jets [17]. In these cases, as explained by previous authors pTl [TBI [T9l [20l [29] . 
the infalling matter first becomes supersonic passing through the outer saddle- 
type sonic point, then the shock forms and it becomes subsonic around a 
center- type sonic point, and finally passing through the inner saddle- type sonic 
point it falls into the black hole. A detailed behavior of such solutions around 
a fast rotating compact object should be investigated for understanding jet 
physics in future. 
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